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In this talk I discuss neutrino oscillations due to possible nonuniversality of the gravitational 



coupling constants. It has been pointed out |jj that a breakdown of the universality of the 
gravitational couplings to different neutrino flavors could lead to neutrino oscillations. In 
• particular the authors of [fj studied the possibility in which the solar neutrinos can be used 
to test this kind of breakdown of the universality. Since the flux of the solar neutrinos is 



relatively small and the energy spectrum is beyond our control, the utility of the solar neutrino 
for this purpose is limited. Here we would like to propose a possible long-baseline experiments 
of neutrino oscillations to test the breakdown of the universality of the gravitational couplings 
to neutrinos [[J H] ||. As we will see, the breakdown of the universality of the gravitational 
kf) 1 couplings to neutrinos of different flavors leads to a violation of Einstein's equivalence principle 
(see, e.g., |J) which states that all the laws of physics must take on their familiar special- 
relativistic forms in any and every local Lorentz frame, anywhere and any time in the universe. 
No consistent theory is known to predict such nonuniversality of the gravitational coupling 
constants, so our motivation is to give an upper bound on such nonuniversality. In the present 
case, it turns out that we can probe the magnitude of the breakdown of Einstein's equivalence 
principle to the order of 10~ 15 , assuming that there are neutrino mixings. Among various 
experiments to test the equivalence principle (see e.g., Ref. for a review), there have been few 
tests of Einstein's equivalence principle for neutrinos ||. The universality of the gravitational 
couplings that I discuss in this talk is of different type from these experiments in the past, so 
our discussions here are complementary to them. 

In this talk we assume that there are two neutrino flavors which have different couplings 
to gravity and that the eigenstates of these different gravitational couplings do not coincide 
with those of the electroweak flavors. Throughout the present discussions we consider neutrino 
oscillations between two flavors for simplicity. Let us suppose that two kinds of neutrinos 
are in a certain background field with the different gravitational couplings, and that they are 
described by the following Lagrangian 
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where we have included mass terms to keep generality, e afl (Gi) (i = 1, 2) are the vierbein fields of 
some background metric with different Newton constants Gi (i = 1, 2), and e(Gi) = dete"(Gj). 
For simplicity we assume that the eigenstates of the gravitational couplings coincide with those 
of the masses. Notice that even if these neutrinos are massless, we cannot rotate these two fields 
so that these are the eigenstates of the electroweak theory, since the gravitational coupling terms 
are not invariant under the rotation in the flavor space. Since the gravitational couplings for 
these two kinds of neutrinos are different, even if we choose a coordinate system in which the 
Dirac equation for V\ in ([!]) becomes the one in a flat space-time, the Dirac equation for v 2 in 
the same coordinate system does not necessarily do so. Thus Einstein's equivalence principle 
is violated in (Jl|). 

The configuration of the long-baseline experiment we will discuss is depicted in Fig. 1, and 
the neutrino beams go underneath the ground. The density of the Earth is not constant, so in 
principle we have to regard the density as a function of the radius |§. When we consider the 
case like the longbase-line experiment in the DUMAND project [TO], where the neutrino beam 
goes deep in the Earth (L =(length of the trajectory of the neutrino beam)~ 6,000 km), this 
is actually the case. However, in the case like the one in the SOUDAN2 project ||10|| , where 
the neutrino beam goes just slightly under the surface of the Earth (L ~ 800 km), we can 
assume that the density is approximately constant. In this talk we will take this assumption 
for simplicity, and the analysis without this approximation is given in [|J. 




Fig. 1 



First let us consider the Dirac equation of left-handed neutrinos without any flavor in the 



interior Schwarzschild background [TT 



{le^D^ - m)il> = 0, 
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where e apL is the vierbein of the interior Schwarzschild metric 
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and is given by 




1 & 



r sin 9. 



(4) 



b(J ab )ip is the covariant derivative acting on a spinor if>, is the spin 



D^ip = (dp - 

connection given by et,uj^ b = d^e^, and a in (f|) is the Schwarzschild radius. 

One characteristic dimensionless parameter in our case is ER, where E is the energy of the 
neutrino, and R is the radius of the Earth. For £7=10 GeV and i?=6,400 Km, ER ~ 3 x 10 23 , 
and derivative terms in the spin connections are all of the order of 1/ER, so we will neglect 
them throughout this talk, a in eq. @ is the Schwarzschild radius of the Earth which is about 
9 mm, so we also expand @ to the first order in a/r. In this approximation the positive energy 
part of the Dirac equation finally becomes 
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where u = Ex is a dimensionless coordinate along the x-axis, t= J R 2 — (L/2) 2 is the distance 
of the trajectory of the neutrino beam from the center of the Earth, and B = El is a very large 
number. 

In case of two kinds of neutrinos described by the Lagrangian ([ID , it is straightforward to 
see that the Dirac equation for ([!]) is given by 
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Here Am 2 



m 



m\ is the difference of the masses, we have defined the Newton potential 



$ = —GMj R on the surface of the Earth, and we have also defined the difference Af — f 2 — fx 
of the dimensionless gravitational couplings of the two neutrino species 
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The equation @ can be easily integrated from u = —EL/2 to u = EL/2. 

Now let us introduce the flavor eigenstates u a , v b of the weak interaction by 

/ cos 9 — sin 9 
I sin^ cos 9 

Then the probability of detecting a different flavor uj, at a distance L after producing one 
neutrino flavor u a is given by 
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This formula applies to the transition between and v T , where no MSW effect |T2] is expected 
to occur. 



In case of the transition between v e and v^, we have to take the MSW effect |12[ into 



consideration, and the Dirac equation is modified as 
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where Gp is the Fermi coupling constant, N e is the density of electrons in the Earth. (|8|) can 
be solved in the same way as before by introducing the variables 



A N (u) cos 26 N 
A N (u) sin 26 N 
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Note that &n does depend on the variable u in this case. It is easy to integrate @, and we 
have the transition probability of detecting v e at a distance L from the source of beams 
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where we have used the fact 6n(u = ^) = 6^{u = — and A at, 9jq are defined through @. 
The integration in the exponent in ([[(]) can be performed numerically 0. 

We have studied the quantity P{y a v b ) for various cases ||. As in the case of neutrino 
oscillations due to masses, we can exclude certain regions in the (sin 2 26, Af) plot. A typical 
figure is given in Fig. 2, which shows the region (sin 2 26, Af) for u^-u r oscillations which can be 
excluded by the proposed longbase-line experiment in the SOUDAN2 project. From Fig. 2 we 
see that we can give an upper limit as small as 10 -14 on the nonuniversality of the gravitational 
coupling constants of neutrinos. In case of the DUMAND project this upper limit could reach 
lO" 15 a. 
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One important feature of this experiment is the dependence of the probability on the detec- 
tion threshold energy of muons, which is, roughly speaking, proportional to the energy of the 
incident neutrinos. To emphasize this aspect, let us make comparison of the energy dependence 
of the probability for various cases. So far we have considered the case where the particle which 
mediates the force between neutrinos is a spin-2 particle, i.e., graviton. In this case the energy 
defined in the Dirac equation 

is given by 



Ej =2 = (l + ^$j= 2 ) \lf + m? (12) 

On the other hand, if the force is mediated by a scalar or a vector particle, then the energy 
would be given by 



E J=0 = Jp ^ + (m + $ J=0 ) 2 



E J=1 = ^F+m? + $ J=1 , (13) 

respectively. Here $j (J = 0, 1, 2) denotes certain potentials for spin-J forces. Examples for 
such forces are a scalar field in which couples only to tau neutrinos and a torsion tensor 
field (i.e., a dual of an axial vector filed) in Jl4j whose eigenstates are different from those of 
the electroweak interaction. It is well-known that the nonrelativistic behavior of scalar forces is 
similar to that of tensor forces, i.e., gravity. In the ultrarelativistic limit, however, the situation 
changes drastically. In this limit we have from fll2f) and ( |13"D 
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where only terms proportional to m and $j could be relevant to neutrino oscillations. So if all 
neutrinos are messless and if neutrino oscillations occur solely due to the presence of $j, then 
the probability of neutrino oscillations would be given by 

Prob(J = 0) ~ sin 2 20 sin 2 (^|pl 
Prob( J = 1) ~ sin 2 26 sin 2 ($ J= iL) 

Prob(J = 2) ~ sin 2 2^ sin 2 (\p\$j= 2 L) . (15) 
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In Fig. 3 we give the energy dependence of the probability in case of gravity. In case of a 
scalar force and a vector force the energy dependence would look like Figs. 4 and 5, respectively. 
From these figures we conclude that the shape of the energy spectrum depends on the spin of 
the particle which mediates the force, and hence it should be possible to identify the spin of the 
particle which mediates a possible new force by looking at the energy dependence. (|i~5l) also 
explains why we have obtained such a severe constraint on the universality of the gravitational 
coupling constants of neutrinos, as the argument of sine in Prob(J=2) in is proportional 
to the energy of the neutrino. We note in passing that ordinary neutrino oscillations due to 
masses is analogous to the case of a scalar force, as far as the energy dependence is concerned. 
In accelerator experiments like those we have proposed in this talk, the length L is much 
smaller and the energy of neutrinos | p\ is much larger than typical quantities in astrophysical 
observations discussed in |7j]. In longbase-line experiments, therefore, while it is easy to test 
gravity, it would be more difficult to detect a new scalar force or a new vector force such as 
those in 



T3j and |3| 
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In this talk we have proposed long-baseline experiments to test the universality of the gravi- 
tational couplings of neutrinos, and we found that we could probe the dimensionless parameter 
Af as small as 1CT 14 or 1CT 15 which is smaller by a few orders of magnitudes than the upper 
limit on a breakdown of the equivalence principle from different types of experiments. Al- 
though we have not evaluated systematic errors in detail, we hope our analysis will stimulate 
and motivate long-baseline experiments in the near future. 
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